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ABSTRACT 

We present a new analysis of a 9-day long XMM-Newton monitoring of the Narrow 
Line Seyfert 1 galaxy Mrk 766. We show that the strong changes in spectral shape 
which occurred during this observation can be interpreted as due to Broad Line 
Region clouds crossing the line of sight to the X-ray source. Within the occultation 
scenario, the spectral and temporal analysis of the eclipses provides precise estimates 
of the geometrical structure, location and physical properties of the absorbing clouds. 
In particular, we show that these clouds have cores with column densities of at least 
a few 10 23 cm -2 and velocities in the plane of the sky of the order of thousands 
km/s. The three different eclipses monitored by XMM-Newton suggest a broad range 
in cloud velocities (by a factor ~4-5). Moreover, two iron absorption lines clearly 
associated with each eclipse suggest the presence of highly ionized gas around the 
obscuring clouds, and an outflow component of the velocity spanning from 3,000 to 
15,000 km/s. 

Key words: Galaxies: AGN — Galaxies: individual (Mrk 766) 



1 INTRODUCTION 

Time variability of X-ray absorption is a common feature 
in Active Galactive Nuclei (AGN). An analysis of column 
density (Nh) variations, performed on a sample of bright 
obscured AGN with multiple hard X-ray (1-10 keV) obser- 
vations (Risaliti et al. 2002) revealed that Nh variability is 
almost ubiquitous among nearby Seyfert galaxies on time 
scales from a few months to a few years. In order to in- 
vestigate shorter variability time scales, two strategies are 
possible: dedicated campaigns of multiple observations, or 
time-resolved spectral analysis of single, long observations. 
In the past few years, both methods have been used to dis- 
cover Nh variations on time scales from a few days down 
to a few hours for a handful of sources: NGC 1365 (Risal- 
iti et al. 2005, 2007, 2009), NGC 4388 (Elvis et al. 2004), 
NGC 4151 (Puccetti et al. 2007), NGC 7582 (Bianchi et 
al. 2009). The physical implications of these results are that 
the X-ray absorber must be made of clouds with linear di- 
mensions of the order of 10 13 -10 14 cm, velocities in excess of 
10 3 km s , and densities n~10 10 -10 11 cm -3 . Assuming that 
these clouds are orbiting around the central black hole with 
Keplerian velocities, the inferred distances from the X-ray 
sources are of the order of thousands of gravitational radii 
(Rg), corresponding to 10 16 -10 17 cm for black hole masses in 
the range from 10 6 to a few 10 7 solar masses. These physical 



parameters are typical of the clouds emitting the broad lines 
observed in the AGN optical and UV spectra. It is therefore 
natural to conclude that the variable X-ray absorber and 
the Broad Line Region (BLR) clouds are one and the same. 

In order to further test this scenario in other sources, 
we started a systematic analysis of all the archival long 
(MOO ks) observations of bright AGN (intrinsic 2-10 keV 
of at least 10 -11 erg s _1 cm~ 2 ) in order to search for oc- 
cultations due to obscuring clouds crossing the line of sight 
to the X-ray source. Our approach is based on a prelimi- 
nary analysis of the hardness-ratio (HR) light curve, with 
the aim of selecting the time intervals with the strongest 
spectral variations, possibly due to occultations, and a sub- 
sequent complete spectral analysis of the selected intervals. 
A detailed description of this method is illustrated in Risal- 
iti et al. 2009A, 2009B for two XMM-Newton observations 
of NGC 1365. 

Here we present a new analysis of a long (~800 ks) 
XMM-Newton observation of the AGN in Mrk 766 
(z=0.0129), which exhibits strong signatures of multiple oc- 
cultations. 

The Narrow Line Seyfert 1 (NLS1) galaxy Mrk 766 is 
highly variable in the hard X-ray (2-10 keV) flux, on time 
scales as short as a few hundred seconds, as typical for 
NLSls (e.g. McHardy et al. 2004). Being one of the X-ray 
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brightest sources of its class, it has been observed inten- 
sively with the main X-ray observatories, such as ROSAT 
(Molendi & Maccacaro 1994), ASCA (Leighly et al. 1996, 
Nandra et al. 1997), BeppoSAX (Matt et al. 2000). More 
recently, higher-quality XMM-Newton observations allowed 
detailed studies of the X-ray variability (Page et al. 2001, 
Boiler et al. 2001, Vaughan & Fabian 2003), and revealed 
the presence of absorption due to highly ionized iron lines 
(Pounds et al. 2003). 

The longest and highest quality X-ray observation avail- 
able for Mrk 766 has been performed with XMM-Newton 
in 2006 over six consecutive orbits, resulting in a total 
elapsed time of about 10 days of quasi-continuous moni- 
toring. A complete time-resolved spectral analysis of these 
data has been presented by Miller et al. (2007) and Turner 
et al. (2007, hereafter T07). In these works, the observation 
is divided in time intervals of 25 ks, and for each interval 
a detailed spectral analysis is performed, testing two dif- 
ferent scenarios: in the first one, the observed spectral and 
flux variability is ascribed to a scattering component; in the 
second one, the dominant changing component is a warm 
absorber partially covering the central source. Both scenar- 
ios can adequately reproduce the observed spectral and flux 
variability. Interestingly, in both models an extra compo- 
nent is needed, consisting of a partial covering by a relatively 
low-ionization cloud with column density in the range 10 22 - 
10 23 cm~ 2 . The best fit values of the column density and 
the covering factor of this component, tabulated in T07, fol- 
low a similar variability pattern in both the scattering and 
absorption scenarios. In particular, in both models these pa- 
rameters have the highest values, and the highest statistical 
significance, during the first orbit and the first part of the 
second orbit, while they are only marginally significant for 
the remainder of the observation. The similar behaviour of 
N//(t) in both scenarios suggests that it is possible to study 
this component in detail independently of the interpretation 
adopted for the underlying spectrum. 

Here we concentrate our analysis on this varying cold 
absorbtion component, trying to isolate its effects from those 
of the other components. We will show that after an appro- 
priate selection of the time intervals the spectral analysis 
strongly suggests the presence of three occultation events, 
due to three different absorbers, with physical parameters 
typical of BLR clouds. 



2 REDUCTION AND DATA ANALYSIS. 

Mrk 766 was observed by XMM-Newton in six consecutive 
orbits from May 23 to June 03, 2005 (the observation in 
the sixth orbit is short (~30 ks) and is not analyzed here). 
The data have been reduced and analyzed following a stan- 
dard procedure, which we already described in other similar 
works (e.g. Risaliti et al. 2009). The only non-completely 
straightforward aspect of the reduction is the choice of the 
rejection level for the background flares. Since the source 
is bright (F(2-10 keV)~2xlO -11 erg cm^s" 1 ), we can use 
a higher threshold than the conservative level suggested in 
the XMM-Newton reduction guide. In particular, since we 
are interested in a time-resolved study, we can accept a rel- 
atively high noise in some time intervals, which otherwise 
would be completely neglected in the analysis. We therefore 



removed only the time intervals where the background flares 
appear to saturate the instruments. This led to the rejection 
of less than 20% of the observing time. Both the temporal 
and spectral analysis were performed using the PN and com- 
bined MOS spectra and light curves. On average, the MOS 
data are more noisy than the PN one, and their contribu- 
tion to the total counts is of the order of 30%. Therefore, 
all the results are dominated by the contribution of the PN 
instrument. 

Our analysis consists in three main steps: (1) a tem- 
poral analysis of the flux and hardness ratio light curves, 
in order to determine the time intervals with significant 
spectral variations; (2) an analysis of the average spectra 
of each of these intervals, in order to determine whether the 
variations are better explained by absorption variability or 
by changes in the underlying spectral shape, and (3) a de- 
tailed time-resolved spectral analysis of the events, under 
the assumption that N_h variability is indeed the driving 
phenomenon. We restrict our analysis to the data at ener- 
gies E>1.5 keV. We note that this choice, while optimal to 
simplify the analysis of the variable component at E>2 keV 
(and therefore, to search for neutral ~Nh variations greater 
than a few 10 22 cm -2 ), prevents us from obtaining a global 
view of the possible reflector/ warm absorber affecting the 
soft part of the spectrum. For an analysis of this compo- 
nent, we refer to T07. 

Each step of the analysis is described below. 

2.1 Temporal analysis 

In Fig. 1 we show the 2-10 keV flux light curve for the 
whole observation, together with the light curve of the (6- 
10 keV) / (2-5 keV) hardness ratio. As expected from pre- 
vious studies (T07), the flux light curve shows strong and 
rapid variability, on time scales shorter than one ks. On av- 
erage, such strong variations are not observed in the hard- 
ness ratio light curve, suggesting that most of the observed 
variability consists of flux changes with a ^constant spec- 
tral shape. The only evident exceptions are present in the 
first and second orbit, where we detect strong changes of the 
hardness ratio, implying variability of the spectral shape. We 
note that the increase of the hardness ratio in these time in- 
tervals is associated to a decrease in the total flux. Therefore, 
these variations may be due either to (1) a hardening of the 
intrinsic continuum, with a pivot point in the hard (>5 keV) 
part of the spectrum, so that the soft flux decreases more 
than the hard one, or (2) to absorption by clouds with col- 
umn density Nh in the interval 10 22 -5xl0 23 cm -2 , which 
would cause strong absorption in the 2-5 keV band, but lit- 
tle variation in the 6-10 keV band. 



2.2 Spectral analysis of the time-averaged 
spectra 

We extracted a spectrum from the third and fourth orbit 
(INT 4 in Fig. 1), where the emission remains constant in 
shape above 2 keV, and determined a model well fitting our 
data, with reduced x 2 ~l> an d no obvious residual features. 
This model consists of a power law continuum, two Gaus- 
sian emission lines at 6.4 keV and 6.7 keV, respectively, and 
a cold reflection continuum, modeled through the PEXRAV 
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Figure 1. Flux (top, in arbitrary units) and hardness ratio HR=F(6-10 keV)/F(2-5 keV) (bottom) light curves of the XMM-Newton 
observation of Mrk 766. For clarity, only EPIC-PN data are shown. 



Table 1. Baseline model in the 1.5-10 keV range - 3 rd and 4*^ 
orbits 
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a : photon index of the primary continuum; b : reflection efficiency, 
defined as the ratio between the normalization of the reflection 
and direct continuum components; c , d , e : peak energy (in keV), 
width (in keV) and equivalent width (in eV) of the iron emission 
lines; * , 3 : 2-10 keV flux and luminosity, in c.g.s. units, adopting a 
(ho, fi./vfj Q;, )=(0.7,0.3,0.7) cosmology. All quoted errors, in this 
table and elsewhere in the paper, are at the 90% confidence level 
for one interesting parameter. 

component (Magdziarz & Zdziarski 1995) in XSPEC. As 
shown in Fig. 2 and Table 1, this simple model is a good 
representation of the data. We note that the reflected con- 
tinuum is quite high (R~2, implying a high covering frac- 
tion by Compton-thick gas) with respect to the equivalent 
width of the narrow iron emission line (EW~35 eV). This 
may be due to either a low iron abundance (a factor of ~2 
under-abundance of iron would be enough to reconcile at a 
90% confidence level the R/EW ratio with that expected by 
George & Fabian (1991) for an inclination of 30 degrees), 



or to an over-simplification of our fit, which reproduces the 
whole curvature with respect to a simple power law above 
~ 3 — 4 keV with this single component, while the actual 
circumnuclear absorber/reflector may be more complex. In 
order to test the consistency of our model, we tried (a) 
fixing the R/EW ratio to the values predicted by George 
Fabian (1991), and (b) replacing the PEXRAV + line com- 
ponents with the self-consistent REFLIONX model (Ross & 
Fabian 2005). In both cases we obtained a fully acceptable 
fit (x 2 =2435/2365 d.o.f., and 2433/2364 d.o.f., respectively). 
A more detailed analysis, such as the one presented in T07, 
can unveil a higher complexity, and the presence of partially 
covering ionized absorbers. However, we stress that our main 
aim here is not a detailed study of these components, but, 
rather, to obtain a reference model of the time-averaged un- 
derlying emission of the source. 

We then added to the analysis three new spectra, obtained 
from the first and second half, respectively, of the first orbit 
(INT 1 and INT 2 in Fig. 1), and from the interval with 
high hardness ratio in the second orbit (INT 3). We fitted 
these spectra with the same model as above, plus an extra 
component consisting of a partial absorber with free column 
density and covering fraction. A visual analysis of the resid- 
uals (Fig. 5) suggested to further add two Gaussian compo- 
nents, in order to fit two absorption features in the 6.7-7.5 
spectral range. The fits were done separately for each inter- 
val, leaving all the spectral parameters free, only requiring 
a constant reflection component through the whole observa- 
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Figure 2. Top two panels: data, best fit model and residuals for 
the spectra obtained from the third (black marks) and fourth (red 
marks) orbit, labeled as INT 4 in Fig. 1. 3 rd panel: residuals with 
respect to the best model of INT 4 of the two spectra extracted 
from the two time intervals INT 1 and INT 2 in the first orbit. 
Bottom panel: same, for the interval INT 3 in Fig. 1. Note the 
difference by a factor of ~10 in the residuals scale among the 
different spectra. For clarity, only data at E>2.5 keV are plotted, 
in order to better show the cutoff-like shape of the residuals. 



tion. The main results are the following: 

(a) In all cases we obtained a good fit (Xred~l)j with the 
partial covering component highly significant in each of the 
three intervals. 

(b) The best fit values of the column densities and covering 
factors of the absorber are different in each interval. 

(c) The best fit slope of the continuum power law is instead 
constant within the errors in all the spectra (including the 
new ones and the reference ones obtained from the third and 
fourth orbit). If we repeat the fits without the partial cov- 
ering absorber (so forcing the continuum slope to adjust to 
reproduce the spectral variations) we obtain a much worse 
fit, rejected at a 99.99% level in all the three fits to INT 1, 
INT 2 and INT 3. This shows that a change in the contin- 
uum slope cannot reproduce the observed spectral variabil- 
ity. Similarly, we also tried to leave the parameters of the 
reflection component free in each fit, and we found that the 
results do not change significantly. 

(d) The two absorption features in the 6.7-7.5 keV interval 
are highly significant (Table 4), and compatible with the 
FeXXV Ka and FeXXVI Ka absorption lines at rest-frame 
energies of 6.70 and 6.96 keV, respectively. We fitted the two 
lines fixing their energies and requiring a common shift. The 
results (Table 4) indicate that the highly ionized absorber 
responsible for these features is in outflow, with velocities 
from ~3,000 to ~15,000 km s _1 . The three outflow velocities 
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Figure 3. Results of the temporal and spectral analysis of the 
first orbit. Upper panel: hardness ratio light curve. This plot is 
just a zoom of the second panel of Fig. 1. Second panel: intrinsic 
2-10 keV flux light curve, normalized to the value measured in 
INT 4. Third panel: column density light curve, obtained from the 
spectral analysis described in the text. Bottom panel: Covering 
factor (CF) light curve. The two pairs of lines (continuous/blue 
and dashed/red) show two possible temporal evolutions of the CF 
of the two individual clouds giving rise to the observed curve (see 
text for details). 



are not compatible with a constant value. Interestingly, the 
absorption features are present only when a cold absorber 
is also present, while they are not detected in INT 4. Since 
the integration time of INT 4 is much longer than in the 
three intervals with evidence of occultations, it may be that 
a variable highly ionized absorber is present also in INT 4, 
but its signatures are averaged out in the total spectrum. 
In order to check this possibility, we divided INT 4 in eight 
intervals, with a fixed duration of 25 ks, and repeated the 
same spectral analysis as for the total INT 4 spectrum. In 
no case we detected an iron absorption line at a 90% sig- 
nificance level. The estimated upper limits for each INT 4 
sub-interval vary between 20 and 30 eV for the single lines, 
well below the level detected during the intervals INT 1-3. 
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Table 2. Spectral analysis of the three high-HR intervals 
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Best fit parameters for the three intervals INT 1, INT 2 and 
INT 3 (Fig. 1). Model parameters: Column density N^f, in units 
of 10 22 cm -2 ; covering factor Cp; slope T of the continuum power 
law; ratio F/F» between the absorption-corrected flux F and that 
measured in the third orbit (INT 4). The continuum reflection 
and the emission lines are fixed to the values shown in Table 1. 



Table 3. Time resolved spectral analysis 
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Model parameters: Column density Nj^, in units of 10 cm" ; 
covering factor Cp, ratio F/F* between the absorption-corrected 
flux F and that measured in the third orbit (INT 4). All the 
other parameters are fixed to the best fit vaues of INT 4, listed 
in Table 1. a : the column density is fixed, due to the complete 
degeneracy with the covering factor when the latter is compatible 
with C^=0. 



We therefore conclude that the ionized absorber responsible 
for the iron absorption lines is physically associated with the 
cold absorber. 



2.3 Time-resolved analysis of the occultations. 

The previous step of the analysis shows that it is possible to 
interpret the spectral variations in the first and second orbit 
as entirely due to occultations by clouds crossing the line of 
sight to the X-ray source. Here we assume that this is the 
correct interpretation, and that no other continuum com- 
ponents are contributing to the spectral variability (while 
we allow for possible intrinsic flux changes). In this way, 
we are able to perform a more detailed spectral analysis on 
shorter time scales, requiring a constant shape of the contin- 
uum components. We divided the first two orbits in several 
short intervals, as shown in Fig. 3 and 4, and performed a 
complete spectral analysis of each interval. 

The length of the intervals is chosen in order to have 




2xl0 4 4x10" 6x10" 

T (sec) 

Figure 4. Same as Fig. 3, for the second orbit of the XMM- 
Newton observation (INT 3 in Fig. 1). 



a balance between the need to sample the changes in hard- 
ness ratio in short time scales, and the need to have enough 
counts in each spectrum to significantly constrain the ab- 
sorber parameters. The continuum spectral parameters 
(spectral slope, reflection component) are those obtained in 
the fit to the third and fourth orbit discussed in the previous 
step of the analysis, and listed in Table 1. The continuum 
flux is free to vary. 

The main results of this analysis are shown in the lower 
panels of Fig. 3 and 4, and in Table 2, and can be summa- 
rized as follows: 

-1. In all cases, we obtain a good spectral fit. This implies 
that the observed spectral variations can be reproduced 
satisfactorily by absorbing clouds moving across the line of 
sight to the X-ray source. 

-2. The absorbing clouds have column densities in the range 
l-3xl0 23 cm~ 2 , while the covering factors vary from zero 
up to ~90% at the peak of the eclipses. 
-3. We interpret the variations observed in the first orbit as 
due to two different clouds, one uncovering the source dur- 
ing intervals 1A to IE, the other covering the source during 
intervals 2A-2E. Since the minimum measured covering fac- 
tor is ~60% (Table 2), the two eclipses overlap in time, in a 
way that cannot be uniquely determined. The bottom panel 
of Fig. 3 shows two possible evolutions of the single occul- 
tations, both resulting in the observed light curve. 
-4. The two iron absorption lines have been included in all 
fits. In principle, this may be a powerful way to investi- 
gate possible variations of the ionized absorber. However, 
the signal-to-noise of the single intervals is not enough for 
a detailed analysis: in all the intervals we find line parame- 
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Table 4. Absorption line parameters 
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Model parameters for the absorption lines, in the three intervals 
discussed in the text. No line is detected in the fourth interval 
(corresponding to the third and fourth orbit, Fig. 1). a : velocity 
in km/s; °: equivalent width in eV of the Fe XXV Ka and K/3 
lines, respectively. 

ters in agreement (within the errors) with the hypothesis of 
no variations apart from the outflow velocity (Table 4). The 
significance of these detections varies considerably among 
the different intervals: in some cases both lines are detected 
at a >4 a statistical significance, in others we can only es- 
timate upper limits. However, this appears to be related to 
the total counts available in the continuum at the line en- 
ergies, rather than to variations in the lines' strength. As a 
consequence, the upper limits in the intervals with no detec- 
tions are of the same order of (or even slightly larger than) 
the values measured in the intervals with a high continuum 
level. For this reason, we will not discuss possible variations 
of the absorption lines further. Let us stress that uncon- 
straining upper limits are found only in some sub-intervals 
of INT 1-3; as we said before, the analysis of INT 1-4 in 
wider time bins is instead very constraining. 

-5. The results listed above do not depend significantly 
on how we fit the variability of the other spectral compo- 
nents. In particular, since our model consists of three main 
components (the direct power law emission, the cold reflec- 
tion, and the iron emission lines), we can assume different 
scenarios: 1) the continuum reflection and the emission lines 
(both the narrow and the broad one) originate from the 
accretion disc, close to the primary X-ray source, and fol- 
low the continuum variations with a time lag much shorter 
than our time bins; 2) the reflection components (continuum 
and narrow line) are emitted far from the primary source, 
and are constant during the whole observation, while the 
broad emission line follows the continuum variations; 3) the 
continuum reflection and both emission lines remain con- 
stant. We repeated our analysis in these three scenarios, 
and found that the results on the continuum occultations 
and on the absorption lines hold with the same significance 
in all cases. Therefore, we do not discuss this issue further. 
We only note, for completeness, that statistically the first 
and second scenario are equally acceptable, while the third 
scenario produces significantly worse fits, in particular with 
strong residuals at the broad iron line energies. 



3 DISCUSSION 

We presented a time-resolved spectral analysis of a long 
(~700 ks) XMM-Newton observation of the NLS1 galaxy 
Mrk 766, suggesting the occurrence of occultations of the 
X-ray source due to circumnuclear clouds. 
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Figure 5. \ 2 residuals for the six intervals considered in Table 3, 
corresponding to the three occultation events (1, 2 and 3A-C), 
the first interval after the third occultation (3D), the remainder 
of the second orbit (3E-F), and the whole third and fourth orbits 
(4). 

While similar events have been seen in obscured AGNs 
(e.g. Wang et al. 2010, Risaliti et al. 2010, Maiolino et 
al. 2010), our analysis demonstrates that such events can 
also be observed also in type 1 sources, which are normally 
unobscured. Indeed, the spectrum of Mrk 766 shows on av- 
erage little obscuration: the eclipses occur in a small fraction 
(~20%) of the observing time, and would be easily missed 
in a global spectral analysis. 

Interestingly, the analysis of T07, where the spectra are 
extracted from fixed 25 ks intervals, provides indications 
of the occultations described here: the N_s values of the 
cold absorption component reported there are on average 
marginally significant, but show local peaks and higher sta- 
tistical significance in the time intervals where we find the 
occultation events. The only difference with our analysis is 
given by the choice of the time intervals, which in our work 
are optimized based on Fig. 1. 

In the following sections we discuss (1) the size and dis- 
tance of the obscuring clouds, as inferred from the observed 
occultations; (2) the internal structure of the clouds, as sug- 
gested by the temporal evolution of the occultations; (3) a 
possible general structure of the circumnuclear medium of 
Mrk 766 accounting for all the results discussed in this pa- 
per. 

3.1 Size and location of the obscuring clouds 

The spectral and temporal parameters of the three occulta- 
tions provide information on the size and distance of the 
clouds from the X-ray source. Assuming that the X-ray 
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source has a linear size of at least 5 gravitational radii, Rg, 
the velocity of the obscuring cloud should be v>2.5xl0 3 
M6 T4 1 \J ACf km s -1 , where M6 is the black hole mass 
in units of 10 6 Mq, T4 is the occultation time in units of 
10 4 s, and ACf is the covering factor variation during the 
eclipse. The black hole mass estimate for Mrk 766, obtained 
from reverberation mapping of H/3, is M6 = 1.8±1.5 (Bentz et 
al. 2009). The durations of the occultations seen in the first 
orbit are of the order of 40 ks, with a covering factor chang- 
ing by ~40%. This implies vi, 2 >800 km s" 1 . We note that 
this value is not precisely estimated, because the end of the 
fist occultation probably overlaps with the beginning of the 
second one. We will discuss this issue in the next Subsection. 
However our conclusions are not affected by this relatively 
small uncertainty. For the third occultation we have instead 
T~ 10 ks and ACf~80%, which gives v 3 >4,000 km s _1 . We 
note that the ratio between the velocities depends only on 
the ratio between the occultation times (and on the assump- 
tion of an axisymmetric cloud shape with respect to the line 
of sight). We have therefore V3~5vi,2 regardless of the un- 
certainties on the size of the X-ray source. The density and 
distance of the clouds from the central black hole can be eas- 
ily derived from the above estimates(assuming Keplerian ve- 
locity). We obtain a density of the order of 10 10 -10 11 cm -3 , 
and a distance of 10 3 -10 4 Rg. 

This simple argument suggests that the X-ray eclipsing 
clouds have the physical properties and location of the BLR 
clouds. 

We note that the above considerations can provide an 
estimate of the size of the X-ray source. The assumption 
made above of a size of 5 Rg can be considered as a "physi- 
cal" lower limit, useful to estimate lower limits for the cloud 
velocities. Since in the above scheme these velocities are lin- 
early dependent on the source size, we can reverse the argu- 
ment, and estimate an upper limit of the size of the X-ray 
source based on the highest physically acceptable cloud ve- 
locity. Considering that the third cloud has an outflow ve- 
locity of ~16,000 km s _1 , we can assume a velocity on the 
plane of the sky of the same order. This implies a linear 
size of the X-ray source of 20 Rg- Adopting instead as a 
"hard upper limit" the highest outflowing velocity detected 
in low-luminosity AGNs (v/c~0.3, Tombesi et al. 2010), we 
would obtain a size of ~100 Rg- Summarizing, the duration 
of the occultations are a direct observational probe of the 
compactness of the X-ray source, with a best estimate below 
~20 Rg, and a maximum size of the order of 100 Rg- 

3.2 Internal structure of the obscuring clouds 

Hints on the internal structure of the eclipsing clouds are 
obtained from the temporal evolution of the model param- 
eters, as listed in Table 2 , 3 and 4. The clearest indication 
that the clouds have a complex structure, with strong ion- 
ization gradients, is provided by the coincidence between 
the iron absorption lines and the cold absorption: while the 
presence of the lines is clear during all the three eclipses, 
no such feature is present in any other time interval, with 
high statistical confidence (Table 4). Moreover, the outflow 
velocities measured in the three occultations are not com- 
patible with each other. This is relevant in two respects: 1) 
it demonstrates that the three occultations are due to three 
different clouds; 2) it reveals a large range in outflow veloc- 



ities among the obscuring clouds, comparable to the range 
spanned by the transverse velocities. We discuss in detail 
these two aspects, and the column density structure of the 
eclipsing clouds. 

1) The three eclipsing clouds. While it was obvi- 
ous from the previous analysis that the third occultation is 
due to a single cloud isolated from the other two, the in- 
terpretation of the occultations in the first orbit was less 
straightforward. In particular, the lowest measured cover- 
ing factor in the first orbit is about 0.6. This means that 
either two different clouds partly overlap at the same time, 
or that a single cloud with a complex structure (a hole, or 
a concave shape) is transiting across the line of sight. The 
association of two different outflowing velocities to the two 
phases of the eclipse seen in the first orbit (INT 1 and INT 2) 
strongly suggests that the former scenario is the correct one. 

If two different clouds are covering the source in the first 
orbit, the slopes observed in Fig. 3 are not directly related 
to the clouds velocites along the plane of the sky. Several 
combinations of the motion of the individual clouds can re- 
produce the observed pattern. In Fig. 3 we show the two 
extreme configurations acceptable in the case of a constant 
slope of the CF time evolution. In both cases, the two clouds 
are simultaneously present along the line of sight in one of 
the two phases. However, we point out that these are just 
two examples of many possible combinations: there is no 
strong physical reason to have a constant slope in the cover- 
ing factor curves (the geometrical shape of the cloud and of 
the X-ray source can easily alter this, as shown, for exam- 
ple, in the case of NGC 1365, Maiolino et al. 2010), and if 
we relax the constraint on the constant slope, there are infi- 
nite ways to reproduce the observed behaviour. In principle, 
this degeneracy may be solved through the observation of 
the absorption lines: the presence of both clouds along the 
line of sight would imply the detection of both line pairs. 
However, this check is not decisive with the available S/N: 
we tried to add a second pair of absorption lines in INT 1, 
with a blueshift velocity of 3000 km s _1 (as measured in 
INT 2), and in INT 2 with a blueshift of 11,000 km s _1 
(as measured in INT 1), and in both cases the additional 
lines are not detected, with upper limits on their equivalent 
widths of the order of 30 eVlj We also repeated this check 
using only intervals IE and 2A (where the contribution of 
both clouds may be stronger) and obtained an upper limit 
of the order of 80 eV, i.e. higher than the typical values in 
Table 4. We conclude that the data suggest that on aver- 
age in both INT 1 and INT 2 one single cloud is dominant 
(at least for the highly ionized component), but we cannot 
exclude a significant overlap in the central part of the first 
orbit. 

The conclusions of this part of the analysis are that (1) 
three separate clouds are responsible of the observed occul- 
tations, each with a different velocity along the plane of the 
sky, and along the line of sight; and that (2) it is not pos- 
sible to uniquely disentangle the path of the two individual 

1 We note that, adding to the complication of the analy- 
sis, the difference between the velocities of the two clouds, 
Av~8,000 km s — 1 , corresponds almost exactly to the difference 
in frequency between FeXXV Ka anf FeXXVI Ka. As a conse- 
quence, a contribution from both clouds to just one of the ob- 
served lines is possible. 
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clouds responsible for the two occultations observed in the 
first orbit. 

2) Ionization structure. The lines detected during 
the occultations are due to hydrogen-like and helium-like 
iron ions. The physical properties (density, column density, 
ionization state) of the gas responsible for these absorptions 
are not easily constrained with the observed lines. The mea- 
sured parameters (equivalent widths and line flux ratios) 
can be obtained with several combinations of column den- 
sity, ionization parameter, and turbulent velocity of the gas. 
Morever, the equivalent widths scale linearly with the (un- 
known) iron abundance. A general treatment of this complex 
problem is presented in Bianchi et al. (2005). Here we only 
stress the two conclusions that can be safely drawn from the 
observed data: the ionization parameter of the gas must be 
quite high (£>1,000 erg cm s _1 ) and the gas column den- 
sity must be higher than a few 10 22 cm -2 . On the other 
hand, the absorption observed in the continuum spectrum 
is due to nearly neutral material. We tried to replace the 
neutral absorption in our models with a ionized one, with 
a free ionization parameter, £ and we obtained again a best 
fit with a neutral absorber, and a limit on the ionization 
£<30 erg cm s _1 . These two phases can be explained as 
due to a "cometary" structure, with a low-density, high ion- 
ization "tail" originating from the high-density cloud head. 
The tail may be either behind the head, with respect to the 
observer, or in front of the head, as in a "real" solar system 
comet (Collin-Souffrin et al. 1988). In the former case, the 
whole cloud, including the dense core, must be in outflow, 
and the tail may be due to friction with a lower density, 
higher temperature gas phase surrounding the cloud. In- 
stead, in the scenario where the tail is in front of the head, 
the dense core can be in a purely Keplerian motion, and 
the outflowing tail may be due to radiation pressure. While 
we cannot distinguish between these two scenarios based on 
the X-ray spectra only, we can safely rule out a spherically 
symmetric cloud structure, without a "head" and a "tail". 
This is due to the requirements of (1) having both the hot 
and cold phases along the line of sight at the same time, 
and (2) a difference between the ionization parameters of 
the two phases of at least a factor of ~40. It is impossible 
to satisfy both these conditions with a spherical cloud: the 
low-density external part, in order to have the required col- 
umn density of at least a few 10 22 cm -2 , must have a linear 
size of at least 10 times that of the cold, dense part. In a 
spherical geometry, this would imply the presence of high 
ionization absorption features for a long time before and af- 
ter the observed eclipses. Regarding the density structure of 
the clouds, we may have either a sudden density drop, or a 
continuous density gradient. In the latter case, intermediate 
ionization states would be present in the cloud. Both these 
scenarios are possible, and our observations cannot distin- 
guish among them, because the tail lies along the same line 
of sight as the neutral head. This hides the effects of possi- 
ble ionization gradients, since the cold absorption removes 
all the marks of a warm absorber, except for the iron-K 
absorption lines, which are the only strong features at en- 
ergies higher than the photoelectric cut off. We note that 
a low-column density tail not covered by the neutral head 
may be present in the third cloud. A low-significance indi- 
cation of such a component may be the marginal detection 
of the absorpion lines in the interval 3D (Fig. 5, Table 3), 




Figure 6. Scheme of the proposed structure of the absorber re- 
sponsible for the observed occultations. Panel 1: structure of an 
eclipsing cloud, as viewed from a direction orthogonal to the line 
of sight. Panel 2: same, as viewed from the line of sight, with the 
ionized tail covered by the neutral head (either behind or in front 
of it). Panel 3: a possible structure of the circumnuclear gas in 
Mrk 766, as viewed from the disc/torus plane. Panel 4: velocity 
components of the eclipsing clouds, as viewed from the disc/torus 
axis. The cloud tail can be either in front (cloud 2) or behind (1) 
the cold core with respect to the source. 



i.e. after the neutral eclipse is over. This line would not be 
statistically significant if detected in an isolated spectrum 
(Ax 2 =6, Table 3) however, given the presence of similar 
lines in the previous spectra, it may be possible that this is 
a real feature. A similar analysis is not possible for the first 
two clouds, given the overlap between the two eclipses, and 
the lack of monitoring during the final phase of the second 
eclipse. 

3) Column density structure of the eclipsing 
clouds. Hints on the column density (Nh) gradients across 
the line of sight can be inferred by comparing the N_» and 
covering factor (Cf) light curves in Fig. 3 and 4. In gen- 
eral, the evidence of Nh changes is not strong (the case of 
constant column density can be rejected at a ~2<r level for 
INT 2, and at a ~3<r level fot INT 1 and INT 3). However, 
we clearly see a general trend of a column density follow- 
ing the covering factor variations: the highest Cf values are 
found together with the highest Nh, and vice- versa. This is 
especially clear in INT 3 (Fig. 4) where we see all the phases 
of the eclipse. We interpret this as an evidence of a column 
denisty gradient across the line of sight, as shown in Fig. 6 
(Panel 2). The possible relation between Cf and Nh may 
be due to the partial degeneracy between these two parame- 
ters in the spectral fits. In order to check this, and in general 
to verify the precision of our error estimates, we calculated 
the Cf-Nh contour plots for several representative intervals 
of the first and second orbit. The results, shown in Fig. 7, 
demonstrate that the degeneracy between the two parame- 
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Figure 7. Covering factor (Cf) - Column density (Net) contour 
plots for several representative intervals. Labels refer to the in- 
tervals as defined in Fig. 3 and 4. 

ters is not severe, and that the single-parameter errors used 
in this work are representative of the real uncertainties. 

At a more speculative level, we notice that the intrinsic 
continuum flux during the third eclipse is anticorrelated with 
Nh and Cf- This may be a simple coincidence, given the 
strong and fast intrinsic variability of this source (Fig. 1), 
or otherwise it could be the indication of the presence of a 
Compton-thick core in the cloud. The effect of this compo- 
nent would be to decrease the observed flux, with no sig- 
nificant change in the spectral shape (except for a small 
flattening of the total spectrum, due to the higher fraction 
of reflected continuum). The possible Compton-thick core is 
also indicated in Fig. 6 with a question mark. 

3.3 A general picture of the circumnuclear 
medium 

The above analysis shows that, within the adopted scenario, 
we can infer several physical and geometrical properties of 
the eclipsing clouds. Here we discuss a general view of the 
circumnuclear region of Mrk 766 accounting for the observed 
occultations. 

The main point that we want to stress is that the ob- 
served clouds are not expected to be "average" Broad Line 



Region clouds. Mrk 766 is a well known Narrow-Line Seyfert 
1, with optical line widths of the order of 1,100 km s _1 
(Grupe et al. 2004) centered at the rest-frame line energies, 
and on average with no absorption along the line of sight. 
This is also the case for the observation discussed here: the 
source remains completely unobscured for ~80% of the time. 
The scenario we propose is described in Fig. 6: the "aver- 
age" broad line clouds form a toroidal structure around the 
central source. The clouds responsible of our occultations 
are out of the main distribution, probably driven by radia- 
tion pressure from the central accretion disc. The proposed 
structure is in agreement with wind-torus models such as in 
Elvis (2000), Proga et al. (2000) or Risaliti & Elvis (2009), 
where a small part of the circumnuclear gas is accelerated 
at outflowing velocities of the order of those observed here, 
in a nearly radial direction, while most of the gas falls down 
into the disk, forming a toroidal stucture. The simulations in 
Proga et al. (2000) show a large range in density (and there- 
fore in ionization) inside the wind region. Since in Mrk 766 
we are not observing a constant wind, it is possible that the 
observed clouds are "outliers" slightly scattered out with 
respect to the main wind stream. 

The clouds detected in Mrk 766 are completely different 
from the ones analyzed in a recent Suzaku observation of 
NGC 1365 (Maiolino et al. 2010), where we were also able 
to observe comet-shaped clouds: in the case of NGC 1365 the 
eclipses are much more common, the motion of the clouds 
is likely dominated by the radial component, and the tails 
lie on the plane of sky. NGC 1365 is a heavily obscured, 
optically type 1.8 source, so we believe that the line of sight 
to the X-ray source is crossed by standard" clouds lying in 
the bulk of the toroidal Broad Line Region. 



4 CONCLUSIONS 

We have presented a time-resolved spectral analysis of 
a long XMM-Newton observation of the NLS1 galaxy 
Mrk 766, suggesting the presence of three occultation 
events due to three different clouds crossing the line of 
sight. Within this scheme, the spectral analysis of the 
single intervals revealed important features of the absorber 
structure: 

1) The obscuring clouds have neautral cores with column 
densities of a few 10 23 cm" 2 , density of 10 10 - 10 11 cm" 3 , 
and linear dimensions of the order of 10 12 -10 13 cm. The 
X-ray source has a similar size, corresponding to a few 
gravitational radii for the estimated black hole mass of 
Mrk 766 of 1.8xl0 6 M . 

2) The detection of Fe XXV and Fe XXVI absorption lines 
during all the three occultations suggests the presence of 
highly ionized cloud tails, outflowing with velocities from 
3,000 to 15,000 km s" 1 . 

Our analysis shows that occultation events, already ob- 
served in several obscured AGN, can happen also in type 1 
objects, which are on average unobscured. When this hap- 
pens, allowing for possible absorption variability is crucial 
both to obtain information on the structure of the X-ray ab- 
sorber (and of the broad emission line clouds) and to perform 
a correct estimate of the parameters of the X-ray spectrum. 

We are currently expanding our time-resolved analysis 
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to several tens of long archival observations of bright AGNs, 
with the aim of increasing the statistics on X-ray occulta- 
tion events, and building a well defined sample in order to 
estimate the frequency of such "black hole eclipses" . 
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